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Abstract 

We point out that formation time of pions produced in heavy ion colli- 
sions modifies the mass spectrum of dileptons produced via 7r + 7r~ — ► e + e~. 
Increasing formation time enhances the production of dileptons with lower 
masses. The effect offers an explanation of a part of the enhanced production 
of dileptons below the p-meson mass as observed by the CERES and HELIOS 
Collaborations at the CERN SPS. 
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I. INTRODUCTION 



Dilepton production in heavy ion collisions is considered as a relatively clear signature 
of the nature of matter produced in heavy-ion collisions. This is one of reasons why the 
enhanced production of dileptons as observed recently in heavy ion collisions by CERES 
[0-0], HELIOS §§ and NA-38 g|7| has received so much of attention. 

An attractive explanation of the peculiar dilepton mass shape is based on the Brown- 
Rho |8| conjecture of decreasing p-meson mass with increasing density and temperature of 
hadronic matter, related to partial restoration of chiral symmetry. 

When included into the transport codes the decreasing p-meson mass together with 
effects due to higher nucleon and A resonancies can lead to a qualitative agreement with 
data U . Good description of dilepton mass spectra in S-Au and S-W interactions has been 
obtained under the same assumption of the behaviour of m p by Li,Ko and Brown [|lO|j . 

The purpose of the present note is to point out that the formation time of pions produced 
in heavy-ion collisions leads also to modifications of the shape of dilepton mass spectrum 
enhancing the production of lighter dileptons with respect to heavier ones. 

The concept of the formation time has appeared a long time ago in connection with 
the problem of how fast is the electromagnetic field accompanying a charged particle recon- 
structed after an abrupt change of the momentum of the particle. Early literature on the 
question can be traced back e.g. from Ref. [|ll|]. In connection with production of secondary 
particles in hadronic collisions the concept has been introduced by Stodolsky [|12| . In analy- 
sis of data [0,0 on multiparticle production in proton - nucleus interactions it has turned 



out [p^ -[T6| that the cascading of secondaries is considerably lower than expected under the 
assumption that a secondary pion is able to interact immediatelly after it has been produced 
in a nucleon - nucleon collision and this has been ascribed to the formation time of secondary 
particles. 

The term "formation time" has somewhat different meaning in different models. In 
what follows we shall understand this term as the time interval between the time of the 
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collision in which the secondary hadron has been produced and the time when this hadron 
is able to interact strongly (with the full cross-section) with other hadrons. For the dilepton 
production in tt + tt~ interaction the formation time is relevant, since the process goes via the 
intermediate p-meson stage 7r + 7r~ — ► p — > e + e~ and the first part of the process is strong. 

We also assume that the formation time of a pion is characterized by t/ which denotes 
the formation time of the pion in its rest frame and in other frames it is Lorentz dilated. 

To see in a very qualitative way the effect of formation time on dilepton spectra consider 
a set of pions produced at t = in a volume of typical dimension L. A pion with four- 
momentum (E,p) is unable to interact before time t = Tf ■ E/m n where 7 = E/m^ is 
the Lorentz time-dilation factor. Within this time the pion traverses the distance d = 
vt = Tf\p\/m n . Faster pions traverse longer distances before interacting. As a consequence 
thereof, interactions of slower pions are prefered and this leads to softening of dilepton 
spectra produced via 7r + 7r~ — > e + e~ . 

This mechanism will be important in situations when each of pions undergoes only a few 
collisions. When the number of collisions, even for fast pions is large, the system thermalizes 
and the spectrum of dileptons corresponds to interactions of fully formed pions in such a 
thermalized system. 



II. SCHEME OF CALCULATION. 

Suppose that the system of pions produced in a heavy ion collision is described by a 
distribution function f(k,x) in phase space, where x = (x,t). The standard expression for 
the rate of dilepton production is (see e.g. fl8fl ): 

' '^Ji^x) f ^%f(k 2 ,x)[v r a e (M 2 ) • 5(M 2 - (h + k 2 f)} (2.1) 



dM 2 d 4 x J (2tt) 3J v x ' ' J (27t) 3 ' 

where a e (M) is the cross-section for 7r + 7r~ — ► e + e~ at given y/s = M, s — {k\ + k^) 2 , 
and v r is the relative velocity. Distribution functions f(k, x) depend on the cross-section a n 
for 7r + 7r~ — > 7r + 7r~ which regulates scattering of pions, and on the value of the formation 
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time. Multiplying and dividing the right-hand side of Eq.(2.1) by a n and integrating over 
d A x we obtain 



R 



<Te(M 2 



dM 2 a w (M 2 ) J (2tt) 3 

which can be rewritten as 



7-1 d 3 k± d 3 k 2 , , 2 , , ; .2, i r2 , 



(2tt) : 



/(A; 2 ,x)[o(M^ - (fci + fc 2 ) 2 ) ■ v^M 2 )] (2.2) 
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where is simply distribution of the number of tctt collisions as a function of the 
mass of the tctt system and R denotes the ratio of the number of tt + tt~ collisions to all tttt 
interactions. The expression is calculated by Monte-Carlo method in a given cascade 
code. The term a^{M 2 ) is the total tttt cross-section including s— and p— waves plus possibly 
higher ones. The term a e (M 2 ) is the 7r + 7r~ — > e + e~ cross-section which can be taken in 
the simplest approximation as 
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This result follows directly from a calculation of the Feynman diagram. Neglecting the 
electron mass m e with respect to yfs = M we obtain 
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The form- factor squared | F n \ 2 can be simply parametrized in the Bright- Wigner approx- 
imation as 



™> 4 P + m 2 p T 2 p 



(2.6) 



1 7r| (M 2 — m 2 ) 2 + m 2 r 2 

In practical calculations we have used a multiple scattering model |19j inspired by that 
used by Humanic |[20|| , see also Ref. pi| . Since tttt interactions are most likely responsible 
for the enhanced production of dileptons with masses below lGeV/c 2 we have studied only 
a simple version of the multiple scattering model in which all secondary hadrons are pions. 
The initial state for the evolution of the cascade is given by set of momenta and positions 



in which pions are created in nucleon - nucleon sub-collisions. This initial state has been 



generated by the Monte-Carlo model of Zavada fllEf . We have used isospin averaged nn 
cross section (see e.g. Appendix in Ref. p2|): 

^r 1 = 7 {l sm26 °° + l smHl + 2 { sm26 > s2e ) ( 2 - 7 ) 

The total cross-section is obtained by integrating over a half of the angular space (due 
to identity of particles) 

ovO) = ~t ( \sin 2 5^ + ^sin 2 5 2 + sin 2 5l) (2.8) 
q \9 9 / 

Here 5f is the phase shift in the partial wave with isospin T and orbital momentum I, 9 



is the CMS scattering angle, s is the CMS energy squared and q = (1/2) ys — Am 2 , denotes 
the CMS momentum. 

Rescattering program follows trajectory of each pion in short time steps and when the 



distance between two pions is smaller than ^a 7T (s)/'K the pions scatter according the cross- 
section given by Eq.(2.7). The scattering angle is determined in CMS frame of two pions 
being scattered and then momenta of pions are transformed back to the global frame of the 



simulation. The scattering is assumed ]H||2"0"|] to take time interval ti of the order of lfm/c, 



which roughly corresponds to h/T p . During this time interval the participating pions are 
unable to interact again. 

In contradistinction to other Monte-Carlo cascade codes, like e.g. [2(| we take initial 



conditions directly from a realistic model |16[ of pion production in nuclear collisions and 



introduce the formation time Tf in the way discussed in Ref. |[24|| . The formation time 
characterizes the formation of pion in its rest frame. 

For times lower than Lorentz dilated formation time (E^/m^Tf a pion created in the 
initial state is unable to interact. The formation time is Lorentz dilated, whereas the time ti 
is fixed in the global frame of simulation which is CMS of Pb+Pb system. Rate of collisions 
per pion varies from 3.7 for Tf = 0.2/m to 1.0 for r = 1.0/m in our simulation. We guess 
that effects of non-locality (see e.g. present also in our rescattering simulation do not 
influence substantially results for dilepton mass spectrum. 



III. RESULTS OF THE SIMULATION 



We have run the rescattering program for three different values of formation time pa- 



rameter Tf using initial state of pion gas generated by cascading generator [16] for Pb+Pb 
160 GeV/n b=7fm collisions. Mass spectrum of dileptons was calculated according to the 
equation (2.3) from the mass spectrum of tttc collisions obtained from the simulation. Cross 
sections given by (2.4) and (2.8) were properly normalized and resulting distribution of dilep- 
tons was rescaled to the 100 MeV bin size, charged multiplicity dN c h/dr] = 300 and yield 
of dileptons per unit of rapidity. In Fig.l we present the shape and magnitude of dilepton 
mass spectrum as a function of the formation time Tf. The interpretation of this dependence 
is simple. Dilepton yields for low masses are almost independent on Tf since slower pions 
are formed rather fast. On the other hand dileptons with masses in the p- meson region are 
suppressed for larger values of the formation time. This is due to the Lorentz dilation factor 
{E^/m^) which permits faster pions to escape from the interaction region with little or no 
rescattering or annihilation. 



CASCUS 1.0E Pb+Pb 160 GeV/n b=7fm 
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Fig.l Dilepton invariant mass distributions obtained by rescattering simula- 
tion for three different values of formation time parameter. 

In Fig. 2 we compare the excess of dileptons as observed by CERES collaboration |26|] with 
our calculations of dilepton production due to nir annihilation as given in Fig.l. Data points 
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and the line which gives the expected dilepton production as extrapolated from proton- 
nucleus interactions are taken from pBfl . Our results presented in Fig.l are added to the 
extrapolated continuous line. 
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Fig. 2 Excess of dileptons from tttt annihilation in hadron gas. Data points 
and background are taken from pf|. 

The comparison shows that a rough and qualitative agreement could be reached for 
formation times Tf ~ 0.5 fm. This comparison (Fig.2) can only serve as a very qualitative 
argument, since the extrapolation procedure usually used to get the background line in Fig.2 
is based on assumption that p-meson dilepton decays and other contributions to dilepton 
production can be simply rescaled from proton-nucleus to nucleus-nucleus collisions. This 
may be quite a rough approximation since p-mesons and other resonancies rescatter in 
Pb-Pb collisions with much larger probability than in proton-nucleon case. We think that 
quantitative analysis can be based only on a detailed Monte-Carlo code containing p-mesons 
and other resonancies in the initial state and taking their rescattering into account. 

IV. COMMENTS AND CONCLUSIONS 

We have shown above that the value of formation time 77 strongly influences the shape 
and magnitude of the mass spectrum of dileptons produced in tttt annihilation. This concerns 
in particular the ratio of low mass dileptons to those in the p-meson region. 
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We have also shown that the shape and order of magnitude of the dilepton production 
due to 7r + 7r~ annihilations roughly corresponds to the excess of dileptons as observed by the 
CERES and HELIOS collaborations at the CERN SPS, both in the low mass and p-meson 
region. Contribution around m ee ~ 500 — 600 MeV/c 2 may still be missing. 

The excess of dileptons as observed at CERN-SPS is of a similar shape as that observed 
by DLS collaboration ^7|. We think that the question whether formation time effects in 



calculations like 28,29 can bring models closer to data deserves a detailed study. 
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